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Signalling pathways are abstractions that help life scientists structure the coordination of cellular 
activity. Cross-talk between pathways accounts for many of the complex behaviours exhibited by 
signalling pathways and is often critical in producing the correct signal-response relationship. For- 
mal models of signalling pathways and cross-talk in particular can aid understanding and drive ex- 
perimentation. We define an approach to modelling based on the concept that a pathway is the 
(synchronising) parallel composition of instances of generic modules (with internal and external la- 
bels). Pathways are then composed by (synchronising) parallel composition and renaming; different 
types of cross-talk result from different combinations of synchronisation and renaming. We define a 
number of generic modules in PRISM and five types of cross-talk: signal flow, substrate availabil- 
ity, receptor function, gene expression and intracellular communication. We show that Continuous 
Stochastic Logic properties can both detect and distinguish the types of cross-talk. The approach is 
illustrated with small examples and an analysis of the cross-talk between the TGF-/3/BMP, WNT and 
MAPK pathways. 



1 Introduction 

Signalling pathway^] are well-known abstractions that help life scientists structure the coordination of 
cellular activity. Interaction between pathways, known as cross-talk, appears to have arisen for several 
reasons, for example to integrate signals, to produce a variety of responses to a signal, to reuse proteins 
between pathways; it accounts for many of the complex signalling behaviours. 

This paper investigates modelling and analysis of pathway cross-talk; a key outcome is the suitability 
of process algebraic operators to modelling cross-talk. We define an approach to modelling based on 
the concept that a pathway is the (synchronising) parallel composition of instances of generic modules 
(with internal and external labels). Pathways are then composed by (synchronising) parallel composition 
and renaming; different types of cross-talk result from different combinations of synchronisation and 
renaming. We use the PRISM modelling language and model checker. The contribution of the paper is 
the following: 

• novel categorisation of types of cross-talk 

• modelling based on pathway module instantiation, internal/external reactions & synchronisation 
over subsets of (possibly renamed) external reactions 

1 we refer simply to pathways henceforth 
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• examples of modelling each type of cross-talk 

• detection and characterisation of cross-talk using Continuous Stochastic Logic 

• application of approach to a case study: cross-talk between TGF-/3/BMP, WNT and MAPK path- 
ways. 

This paper is organised as follows. The following section outlines the background to signalling 
pathway cross-talk and defines a novel categorisation based on examples in the literature. In section [3] 
we describe the modelling approach. In sectionHwe show how each type of cross-talk is modelled for an 
example pathway and how we can detect and characterise the cross-talk using CSL. Section [5] describes 
the case study, cross-talk between TGF-/3/BMP, WNT and MAPK pathways, and gives some results. 
Section [6] contains a discussion of our overall approach and Section previews related work. Conclusions 
and directions for future work are in section 

Throughout the paper the following notation is used. Transformation (e.g. protein X turns into 
protein Y) is denoted by a solid line with an arrow. Catalysis (increase in the rate of a reaction) is 
denoted by a dashed line with an arrow. Inhibition (decrease in the rate of a reaction) is denoted by a 
solid line with a blunt end. Finally, we distinguish between inactive proteins and active proteins rather 
than the various mechanisms by which a protein changes state. An active protein is decorated with *, e.g 
X is inactive and X* is active. This notation is illustrated in Figure (Ha). 



3 ^ Transformation > ^) 



Catalysis > 

Inhibition 1 

Protein OD 
Activated Protein 



Figure 1: a) the notation used throughout this paper for arcs and nodes, b) an example of a 2-stage 
signalling cascade in which the activated protein X catalyses the activation of protein Y. 



2 Signalling Pathways & Cross-talk 

Pathways are the mechanism by which a cell receives a signal and produces the appropriate cellular 
response. A cell can respond to many different signals, and the same signalling pathway can produce 
different responses based on different signals. For example, the well-studied MAPK/ERK pathway can 
respond with cellular proliferation or cellular differentiation depending on the type of growth-factor 
present (T). Many of the reactions involved in signalling pathways are enzyme catalysed protein activa- 
tions, often arranged in a "signalling cascade". In such a cascade, the activated protein on one "level" is 
the enzyme for the activating reaction of the next "level", as shown in Figure QJb). 

Signalling pathways were first thought to be a linear series of reactions, but more recently, detailed 
understanding of these pathways shows that they are non-linear Q. The series of reactions forming 
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the pathway can, for example, diverge or interact upstream/downstream in the chain of reactions form- 
ing a feedback/feedforward loop. Laboratory techniques measuring the concentration of proteins or 
RNAs/DNAs in the cell under different conditions now enable biologists to build appropriate abstrac- 
tions of pathways. This is an incremental process further complicated by the lack of clear definition 
of what constitutes a single pathway. It is an inexact science to define the boundaries between path- 
ways. Often, the boundaries are simply drawn such that the abstraction of the pathway can explain the 
biological data. 

The term cross-talk was first applied to electronic circuits to describe a signal in one circuit having 
an undesired effect on another circuit Q. Cross-talk in this setting is a design flaw: the electronic circuit 
has been specified and built, and has resulted in an undesired interaction between signals, called "signal 
interference". Biochemical cross-talk J4) is an interaction between signals flowing through two or more 
signalling pathways in a cell, however, this is not necessarily indicative of signal interference. 

2.1 Types of Cross-talk 

Cross-talk can occur at all stages of signal propagation through a pathway. Although there is some dis- 
cussion of types of cross-talk Q, there appears to be no universal categorisation in the literature. Here, 
we define five categories of cross-talk: signal flow, substrate availability, receptor function, gene expres- 
sion and intracellular communication. The five categories are illustrated in Figure [2] and are discussed in 
more detail below, with reference to an indicative example. We note that four of the five categories are 
alluded to in O but are not made specific. 

2.1.1 Signal Flow Cross-talk 

Signal flow cross-talk between two pathway occurs when a molecular species in one pathway affects 
the signal flow (rate of protein activation) in another pathway, shown in Figure (2a). The species affects 
the signal flow by altering the rate(s) of the activation or deactivation reactions in the other pathway. 
Typically, this interaction occurs in the cytoplasm and affects the rate of the downstream signal. For 
example, in JU there is signal flow cross-talk between the MAPK and Integrin signalling pathways. 
Activation of the Integrin pathway enhances signalling through the MAPK pathway by increased rate of 
activation of key proteins in the pathway. 

2.1.2 Substrate Availability Cross-talk 

Substrate availability cross-talk occurs when two pathways compete for one or more common or homolo- 
gous proteins (proteins that perform the same function), shown in Figure0b). For example describes 
two pathways that compete for activation of the MAPK cascade. The pathways share the MAPKKK 
protein STE1 1 and have homologous MAPKK and MAPK proteins. 

2.1.3 Receptor Function Cross-talk 

Cross-talk can occur at the level of the pathway receptor, shown in Figure (2c). The receptor's ability 
to detect a ligand can be affected by cross-talk, for example, being inhibited (thus slowing or blocking 
signal propagation) or activated (thus producing signal propagating in absence of the ligand). In (H other 
signalling pathways can activate the estrogen receptor in the absence of the estrogen ligand. 
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2.1.4 Gene Expression Cross-talk 

Cross-talk can occur at the level of gene expression within the nucleus, shown in Figure 0d). In this 
case, there is an interaction between pathways in terms of genes that will be expressed or repressed, for 
example, through activating/deactivating transcription factors (TFs) or increasing/reducing the number 
of TFs. In (9) two pathways contain cross-talk within the nucleus. One pathway contains a transcription 
factor called GR that resides outside the nucleus in its inactive state. Upon activation (signalling), GR 
relocates to the nucleus and represses the transcription factor NF-/cB that is activated through another 
pathway. 

2.1.5 Intracellular Communication Cross-talk 

Finally, signalling pathways can cross-talk using the less direct manner of intracellular communication, 
shown in Figure (2e). Instead of physical interaction between the proteins that comprise the pathways, 
one pathway can release a ligand that activates another pathway. In Q the TGF-/3/BMP and WNT 
pathways reciprocally regulate the production of their ligands. There is some contention in the literature 
as to whether this is genuine cross-talk: the interaction is less-direct than other types of cross-talk and 
involves lengthy processes such as gene expression and ligand excretion. 



3 Modelling 

Our aim is to build models of pathways and their cross-talk in a modular fashion, whilst demonstrating 
our cross-talk categorisation. 

The literature contains many examples of modelling biological systems based on an unstructured set 
of equations, for example Ordinary Differential Equations (ODEs) JTUl . However, flat equations have 
disadvantages such as lack of structure and modularity - there is no accepted way of dealing the with 
composition of these models. Several formal languages from Computer Science have well-understood 
notions of structure. 

We choose a state-based formalism, the PRISM modelling language ifTTTl . with modules, renaming 
and synchronous communication between modules. The semantics is given by continuous-time Markov 
chains (CTMCs), which can be analysed using the PRISM model checker with Continuous Stochastic 
Logic (CSL). We give a brief overview of the language. 

The language is based on "reactive modules", each of which can contain local variables. The vari- 
ables are updated by the execution of commands that have the following syntax, 
[label] guard -> rate : update. statement. Commands are only executable when their guard 
becomes true, and any label synchronises as required. The rate is used to build the underlying CTMC, 
providing both the probability and timing information of the state after update. statement executes. 

Modules can be composed concurrently, synchronising (multiway) on the commands whose labels 
occur in the synchronisation set, L. For example, given modules Ml and M2, and set of labels L, 
Ml |[L]| M2 denotes the concurrent composition of Ml and M2, synchronising on all labels in L. If 
the label set is omitted, Ml 1 1 M2, then Ml synchronises with M2 on the intersection of labels occurring 
in Ml and M2. PRISM also allows renaming of labels, denoted thus Ml {oldJabel ^— newJabel}, and 
hiding, denoted thus M / {labell, . . . , labeln}. Hidden labels are not available for synchronisation. 
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a) Signal Flow 
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Figure 2: An example of each of the five types of cross-talk: a) a pathway up-regulates signal flow 
through another pathway, b) two pathways compete for a protein, c) a pathway activates the receptor of 
another pathway in the absence of a ligand, d) two pathways have conflicting transcriptional responses, 
e) a pathway releases a ligand for another pathway. The following notation is used: L1/L2 - ligands, 
R1/R2 - receptors, W/X/Y/Z - Proteins, TF - Transcription Factor. 



3.1 Generic Modules 

We define a pathway module to be a behavioural pattern within a pathway. For example, commonly 
occurring pathway modules are: receptor, 3-stage cascade and gene expression (Figure [3]). 

We represent these by generic modules in PRISM as follows. We adopt the reagent-centric modelling 
style, as first presented in lfT2l . Local variables contain the concentration of proteins in a particular 
state (e.g. active, inactive). Labelled commands change protein concentration according to biochemical 
reactions. Reactions are considered to be external or internal. The former denote behaviour that can 
coordinate (or be coordinated with) behaviour in other modules, they are available for synchronisation. 
The latter reactions are hidden and are not available for synchronisation. Concentration is modelled by 
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Figure 3: Three generic pathway modules: receptor, 3-stage cascade and gene expression. There are 14 
external reactions (el,... ,el4) denoted by black arcs, 2 internal reactions (/l,/2) denoted by gray arcs 
and 6 species that are present in the initial state denoted by shaded ellipses. 



discrete, abstract levels, as defined in (13 1. Here, we consider strong abstractions, usually two levels 
{0,1}, and unit reaction rates. Due to the ligand production reaction e4, ligand has three levels {0,1,2}, 
and the rate of ligand-receptor binding reflects the level of the ligand (law of mass-action). 
The PRISM modules for receptor, 3-stage cascade and gene expression are the following. 

module Receptor 

R : [0..1] init 1; LI : [0..2] init 1; RIActive : [0..1] init 0; 

[il_l] Rl = 1 & LI >= 1 & RIActive = -> L1:(R1> = 0) & (LI' = 0) & (RIActive* = 1); 

[el_l] RIActive = 1 -> 1: (RIActive ' = RIActive); 

[e2_l] RIActive = 1 -> 1: (RIActive* = RIActive); 

[e3_l] Rl = 1 & RIActive = -> 1 : (Rl 5 = 0) & (RIActive ' = 1); 

[e4_l] LI < 2 -> 1:(L1' = LI + 1) ; 
endmodule 



module Cascade3 

Xllnactive : [0..1] init 1; XIActive : [0..1] init 0; Yllnactive : [0..1] init 1; 
YlActive : [0..1] init 0; Zllnactive : [0..1] init 1; ZlActive : [0..1] init 0; 

[e5_l] Xllnactive = 1 & XIActive = -> 1 : (Xllnactive ' = 0) & (XIActive 5 = 1); 
[e6_l] Yllnactive = 1 & YlActive = & XIActive = 1 -> 1 : (Yllnactive 5 = 0) & (YlActive* = 1) ; 
[e7_l] Yllnactive = 1 & YlActive = -> 1 : (Yllnactive > = 0) & (YlActive* = 1); 
[i2_l] Zllnactive = 1 & ZlActive = & YlActive = 1 -> 1 : (Zllnactive * = 0) & (ZlActive* = 1) ; 
[e8_l] XIActive = 1 -> 1: (XIActive* = XIActive); 
[e9_l] XIActive = -> 1: (XIActive* = 1); 
[el0_l] ZlActive = -> 1: (ZlActive* = ZlActive); 
[ell_l] ZlActive = 1 -> 1: (ZlActive' = ZlActive); 
[el2_l] Xllnactive = 1 -> 1 : (Xllnactive ' = 0); 
endmodule 



module GeneExpression 

Genel : [0..1] init 1; Proteinl : [0..1] init 0; 



[el3_l] Genel = 1 & Proteinl = -> 1: (Genel* = 0) & (Proteinl' = 1) ; 
[el4_l] Proteinl = 1 -> 1: (Proteinl' = 0); 
endmodule 
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We will treat these modules as generic, that is, we instantiate them (strictly, duplicate and rename in 
PRISM) for multiple occurrences. For brevity, we adopt the following convention. For generic module 
M, Mi denotes an instance of M with every variable and label renamed by an indexed form. For example, 
variable v becomes v\ in module M\ . 

3.2 Pathways 

A pathway is a parallel composition of instances of the generic modules, renaming labels to coordinate 
synchronisation within the pathway. 

Definition 1. Let G be a set of generic modules. A pathway P has the form f\X\ \[L\]\ ... | [L n -i] I fnX n 
where X\ ... X n are instances of modules in G, f\ . . . f n are compositions of renaming s and hidings and 
L\ ... L n _\ are labels. 

As an example, consider the expression of a pathway comprising the active receptor catalysing the 
activation of protein X and active protein Z catalysing the expression of Gene. This is defined by: 

P\ — Receptor\ / {il\} {e\\ 4— e5i} |[e5i]| Cascadei\ / {i2\} {ell\ ^— e\3\} 
|[el3i]| GeneExpression\ 

The Receptor i and Cascade?* \ modules synchronise on e5\, and Cascade? \ and GeneExpression\ 
synchronise on e\?\ ; in both cases these are the only external labels that occur in both modules. Internal 
reactions i\\ and i2\ are hidden using the / operator. 

3.3 Composition of Pathways 

In a similar way, pathways can be composed, synchronising on external labels. 

Definition 2. Given two pathways P\ and P2, with sets of external labels ext{P\) and extiP^) resp., if 
ext{P\) n ext{Pj) — {} then the pathways are independent, otherwise there is crosstalk. 

As an example of independent pathways, consider P\ and P2 (shown in Figure 0)), where 
P2 = Receptor^ / {^2} {^2 <— £$2} \ [ e ^2] \ Cascade^ I {^2} {^1 12 ^~ ^^2} 

I [el 32] I GeneExpressioni 
We synchronise Pi and P2 over the unused external reactions, i.e. the reactions that are otherwise 
never involved in synchronisation. In this case, P\ \ \U]\Pi where U = {e2\ , e?\ , eA\ , el\ , e%\ , e9\ , elOi , 
<?12i, <?14i, e22, e?2, ^2, el 2, ^82, ^2, elCb, <?122, el42}. Since the labels in U occur in only one 
pathway they cannot synchronise and their corresponding actions never execute. 

P\ and P2 are independent pathways, since ext(P\) fl extfa) = {}. We can also demonstrate in- 
dependence through CSL properties. For example, the following CSL property expresses that it is not 
possible to activate P\ without activating receptor R\ or activate P2 without activating receptor R2. 
p x \[U}\ P 2 \= P< [F{RActivei = 0A Protein x = 1) V F(RActive 2 = 0A Protein 2 = 1) ] 
In the next section we consider the case where the intersection is not empty; specifically, we consider 
examples of the five possible types of cross-talk between P\ and P2. 

4 Cross-talk models 

We now consider the five possible types of cross-talk between Pi and P2. In the previous section 
we synchronised P\ and P2 over U, the set of unused reactions. We now synchronise over E and U, 
P\ \[E, U} \ P2, where E = ext(P\) fl ext(P2), i.e. the common external labels between P\ and P2, which 
may involve renaming. For each type we give a general description and then an explicit example. 
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Figure 4: The two pathways P\ and P2 each comprising three instances of the generic pathway modules 
receptor, 3-stage cascade and gene expression. External reactions are denoted by black arcs, internal 
reactions by gray arcs and species that are present in the initial state by shaded ellipses. 

We say that there is competition for a protein if the both pathways modify the protein such that they 
become unavailable to the other pathway. An enzyme that is shared between two pathways could be 
considered competition, however the duration of the enzyme-substrate complex is extremely short. As 
such, we ignore this step in our approach. 

Signal Flow Cross-talk P\ \[E, U] \ P2 has signal flow cross-talk if E only contains labels from two 
cascade modules or labels from a cascade module and catalysis/inhibition labels from a receptor module, 
in either case not setting up competition for a protein between the pathways. Example: provide an 
additional route to the activation of Y\ through the X 2 * enzyme. Synchronise ^82, the enzymatic activity 
of X 2 *, with el 1 the alternative route to activate Y\, by renaming e8 2 to el\ and synchronising between 
pathways onel\. 

P X \[el u U]\ P 2 {e8 2 <- eli} 

where U = {e2\, e3\, e4\, eS\, e9\, elOi, e\2\, el4\, e2 2 , e3 2 , e4 2 , el 2 , e9 2 , el0 2 , el2 2 , el4 2 }. 

Substrate Availability Cross-talk Pi \[E, U]\P 2 has substrate availability cross-talk if E only contains 
labels from cascade modules or labels from a cascade module and catalysis/inhibition labels from a 
receptor module, in either case setting up competition for a protein between the pathways. Example: 
make both pathways compete for the activation of protein X\ catalysed by their respective receptors PI* 
and P2*. Within P 2 , synchronise e2 2 , the enzymatic activity of R 2 *, with e9 2 , the production of X 2 *. 
Between pathways, synchronise e\2\, the degradation of X\, with the e9 2 . This produces a new reaction, 
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X\ -> X 2 * with R 2 * as the enzyme. Also, synchronise between pathways on e5 2 to block the following 
reaction, X 2 ^> X 2 * with P2* as the enzyme. 
Pi {ell! <- e9 2 } \[e9 2 , U]\ P' 2 

where U = {e2\, e3\, e4\, el\, <?8i, e9\, elOi, el4\, e3 2 , e4 2 , e5 2 , el 2 , <?8 2 , el0 2 , el2 2 , el4 2 } 
and/2 — Receptor 2 / {il 2 } {el 2 e5 2 , el 2 ^— e9 2 } \[e5 2 , e9 2 ,]\ Cascade3 2 / {i2 2 } 
{ell 2 <(— ^132} I [el3 2 ] I GeneExpression 2 . 

Receptor Function Cross-talk Pi U]\ P 2 has receptor function cross-talk if E only contains 
labels from a receptor module and catalysis/inhibition labels from either a receptor module or a cascade 
module. Example: provide an alternative route to activate receptor R 2 by the enzyme Xi*. Synchronise 
e3 2 , the alternative route to active R 2 , with eS\, the enzymatic activity of X\*. 
Pi {e% x <- e3 2 } \[e3 2 , U]\ P 2 

where U = {e2i, e3\, e4\, el\, e9\, elOi, e\2\, e\4\, <?2 2 , e4 2 , el 2 , <?8 2 , e9 2 , el0 2 , e\2 2 , el4 2 }. 

Gene Expression Cross-talk P\ \[E, U] \ P 2 has gene expression cross-talk if E only contains labels 
from a gene expression module and catalysis/inhibition labels from a cascade module. Example: connect 
the inhibiting activity of Z2* to the expression of Gene\. Synchronise e\0 2 , the inhibiting activity of Z2*, 
with e\3\, the expression of Gene\. 

P\ \[el3 u U]\P 2 {el0 2 <- e\3 x } 

where U 

{e2\, e3\, e4\, el\, e8i, e9\, elOi, e\2\, e\4\, e2 2 , e3 2l e4 2l el 2 , e% 2 , e9 2 , e\2 2 , el4 2 }. 

Intracellular Communication Cross-talk P\ \[E, U] \ P 2 has intracellular communication cross-talk 
if E only contains a protein degradation label from a gene expression module and a ligand production 
label from a receptor module. Example connect the degradation of Protein\ to the production of L 2 . 
Synchronise e\4\, the degradation of Proteinic with e4 2 , the production of L 2 . 
Pi {el4 1 <- e4 2 } \[e4 2 , U}\ P 2 

where U = {e2\, e3\, e4\, el\, e%\, e9\, elOi, <?12i, e2 2 , e3 2l el 2l e% 2 , e9 2l el0 2l el2 2l el4 2 }. 
4.1 Detecting Cross-talk 

We perform Continuous Stochastic Logic (CSL) model checking in PRISM to detect the presence of 
cross-talk. For each of the 5 cross-talk models we compute the probability of three properties and note 
that cross-talk is detected by a change in probability compared with the independent model. The three 
properties are as follows. 

Competitive Signal Flow (Pi): probability of signal flow through Pi before P2 

P =? [F{Protein\ = 1 A Protein 2 = 0) ] 

Time-dependent Signal Flow (Pi): probability of signal flow through Pi within 3 time units 

P=? [F< 3 (Pwteini = 1) ] 

Time-dependent Signal Flow (P 2 ): probability of signal flow through P2 within 3 time units 

P=? [F< 3 (Protein 2 = 1) ] 
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The competitive signal flow property detects the presence of 3 of the 5 types of cross-talk. The 
probability of this property is significantly greater with signal flow cross-talk, due to the extra signal 
flow from P2, and significantly less with gene expression cross-talk, due to inhibition of gene expression 
by P2. There is a decrease in the probability of this property with receptor function cross-talk, however 
the decrease is small due to the likelihood that the cross-talk is initiated after P2's receptor becomes 
activated. The competitive signal flow property does not detect the substrate availability cross-talk, as 
the cross-talk has an identical effect on each pathway. Intracellular communication cross-talk is also not 
detected, as the cross-talk occurs after protein\ is expressed (the property does not refer to behaviour 
past this point). 

To detect substrate availability and intracellular communication cross-talk we need to check the prob- 
ability of independent signal flow through a pathway (hence, not in relation to another pathway). This 
is accomplished by checking the probability of signal flow through a pathway within a time bound, in 
this case (arbitrarily) chosen to be 3 time units. Substrate availability cross-talk has the effect of equally 
decreasing the time-dependent signal flow through both Pi and P2 due to the competition for a limited 
protein. Intracellular communication cross-talk has the effect of increasing only the signal flow through 
P2, however the effect is marginal (5th decimal place). It is interesting to note that in Section 12. 1.51 we 
identified intracellular communication cross-talk as a source of contention amongst the community and 
it proves difficult to detect in our analysis. 

4.2 Characterising Cross-talk 

We now define 5 CSL properties, each of which characterises a type of cross-talk and thus holds only in 
the respective model defined previously. The properties are simple liveness or safety properties and do 
not exploit the rate information in the model. These properties could equally be written in Computational 
Tree Logic, replacing the probabilistic operator with the universal (A) and existential (E) operators as 
appropriate. 

Signal Flow Cross-talk (flow from P2 to Pi) it is possible to activate Pi without activating receptor Pi 

P >0 [F(RActivei = A Protein x = 1) ] 
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Substrate Availability Cross-talk (Pi and P2 compete for a protein) it is not possible to activate 
both P\ and P2 (i.e. the pathways compete for a limited protein) 

P< [F(Proteini = 1 A Protein 2 = 1)] 

Receptor Function Cross-talk (Pi activates /Vs receptor) it is possible to activate the receptor R2 
without using the ligand L2 

P >0 [F(RActive 2 = 1 A L 2 = 1) ] 

Gene Expression Cross-talk (P2 inhibits Pi 's gene expression) it is not possible to express protein\ 
if the signal has already passed through P2 

P< [F{Protein\ = 1) {Protein\ = A Protein 2 = 1}] 

Intracellular Communication Cross-talk (Pi expresses /Vs ligand) it is possible to use and replen- 
ish ligand L2 

P >0 [ (L 2 = 1) A (L 2 = 1) U ( (L 2 = 0) A (L 2 = 0) U (L 2 = 1) ) ] 

5 Case Study: TGF-/3/BMP, WNT and MAPK pathways 

We applied our approach to a prominent biological case study of the cross-talk between the TGF-/3/BMP, 
WNT and MAPK pathways. Details are taken from [5 ] and from discussions with a domain expert lfT4ll . 
We consider the behaviour of the independent pathways and the various types of cross-talk. We note that 
the effects of cross-talk are not discussed in Q. 

Our model of the pathways and their cross-talk is shown in Figure To apply our module approach 
we need to expand our set of modules to: receptor, protein activation, 2-stage cascade, 3-stage cascade, 
translocation, protein binding and gene expression. This is a natural extension of our approach, and 
the extra modules act in a similar manner to the modules that have been discussed. Using the formal 
descriptions of cross-talk (Section 0)), we have identified three types of cross-talk in this model. 

We measure the output of the TGF-j8/BMP pathway by the activity of the expression of Genes to 
Proteins. We use the following properties to compare the effects of cross-talk: \\f\, the eventual expression 
of Genes, and 1//2, the time-dependent expression of Genes (within 5 time units). 
y/i = P =? [F (Proteins =1)], Vi = P=? [ F< 5 (Proteins =1)] 

Independent Pathways With independence, the activation of the TGF-j8/BMP pathway leads to gene 
expression within 5 time units, y/ 2 , with probability 0.47 and eventual gene expression, y/i, with proba- 
bility < 1 due to the inactivation of the receptor. 

TGF-/3/BMP and MAPK Cross-talk There are two types of cross-talk between the TGF-/3/BMP and 
MAPK pathways. Signal flow: MAPK* proteins slow signal flow through the TGF-j6/BMP pathway 
by deactivating the R-Smads and degrading SmadA. Gene expression: the TF* and AKT* proteins 
upregulate gene expression in the TGF-j8/BMP pathway. Note that the appearance of the AKT and 
PI3K proteins in the MAPK pathway indicates an implicit cross-talk with the AKT and PI3K pathways 
respectively. The inclusion of cross-talk with the MAPK pathway can both provide alternative gene 
expression routes and block the TGF-jS/BMP route, overall causing the probability of gene expression 
within 5 time units, y/2, to marginally decrease to 0.73. The probability of eventual gene expression, 
is 1 due to the consistent routes through the MAPK pathway. 
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MAPK TGF-B/BMP WNT 




Figure 5: Cross-talk between the TGF-jS/BMP, WNT and MAPK pathways. Species that are present in 
the initial state are denoted by shaded ellipses. 

TGF-j8/BMP and WNT Cross-talk There are three types of cross-talk between the TGF-/3/BMP and 
WNT pathways. Signal flow: the Srnadl* protein degrades fi-Catenin and the axin protein degrades 
SmadiQ Gene expression: the fi-Catenin protein upregulates gene expression in the TGF-/3/BMP path- 
way. Intracellular communication: the WNT pathway can cause the production of a ligand for the 
TGF-/3/BMP pathway, and vice-versa. The inclusion of cross-talk with the WNT pathway can both pro- 
vide an alternative route to gene expression and inhibit Smadl which can inactivate the receptor for the 
TGF-/3/BMP pathway. Overall this causes the probability of gene expression within 5 time units, y/2, 
to marginally increase to 0.76. The probability of eventual gene expression, y/i, is still < 1 due to the 
degradation of the fi-Catenin protein. 



TGF-/3/BMP, WNT and MAPK Cross-talk The cross-talk between all three pathways is the union 
of the two cross-talk scenarios above. The effect of both WNT and MAPK cross-talk to the TGF-/3/BMP 
pathway is additive. The probability of y/2 has risen to 0.88, compared with the single cross-talks of 
WNT and MAPK with probability 0.76 and 0.73 respectively. The inclusion of the MAPK cross-talk 
provides consistent routes to gene expression and hence the probability of \\f\ is 1. 



2 this cross-talk is discussed further in Section [6l 
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6 Discussion 

Categorisation The categorisation of signal flow cross-talk between the WNT and TGF-/3/BMP path- 
ways in the previous section may be open to discussion. In Figure \5\ Axin degrades Smadl. Rather than 
being part of the signal flow from the TGF-j8/BMP receptor to gene expression, Smadl actually deacti- 
vates the TGF-/3/BMP receptor. One could consider this cross-talk as receptor function cross-talk with 
an intermediate {Smadl). Our approach categorises this as signal flow because the set of labels involved 
is an catalysis label from a protein module and a degradation label from a protein activation module. 

Quantitative Detail We have demonstrated our approach in this paper on models with a low level 
of quantitative detail. As such, the probability values resulting from CSL model checking can only 
be used to compare between models. However, with more quantitative detail, further interpretation of 
our analysis results would be possible. For example, the properties concerning the probability of time- 
dependent gene expression between cross-talk models would become a meaningful assessment of the 
strength of the cross-talk. 

Feature Interaction There may be an interesting analogy with feature interactions in telecommunica- 
tions and software systems. Features, or services, in these systems are additional functionality (additional 
to the core). They are often added incrementally, by various developers, at various times (e.g. due to 
deregulation). A possible consequence is interactions between the new features themselves, or with the 
core system, causing some features or the core to behave in new, sometimes undesirable ways. An open 
question is whether techniques developed to model and detect features and interactions may be applicable 
to pathway cross-talk. Moreover, a common problem is lack of universal definition of pathway/feature; 
it would be interesting to investigate if concepts such as the feature construct of lfT5l would be useful 
in the pathway paradigm. We note that many approaches to interaction detection are based on temporal 
logic descriptions of behaviour. 

7 Related Work 

The literature contains a limited number of applications of computational techniques to the study of 
signalling pathway cross-talk. 

There are several examples whereby a proposed method of cross-talk between two pathways is ex- 
pressed and analysed as a computational model. For example (T6) and ifTTll use Ordinary Differential 
Equations (ODEs) to model the cross-talk between the MAPK pathway with the AKT and PKC path- 
ways respectively. The work in lITSl is similar, however the modelling technique used to model the 
cross-talk between the growth factor pathway and the Integrin pathway is stochastic networks. Further- 
more, a more formal notation of Petri nets has been used in llT9l to model the cross-talk between the 
pathways involved in Apoptosis decision-making. 

The computational analysis of cross-talk models have produced some interesting results in the lit- 
erature. For example, Q has produced a model containing substrate availability cross-talk between the 
hyperosmolar and the pheromone MAPK pathways. The question the authors try to answer is how these 
pathways maintain signal specificity given that they share common proteins. The authors analyse two 
models which can account for the signal specificity, one which contains mutual inhibition between path- 
ways to limit signal bleed-through and one which contains scaffold proteins. In ifTTIl the computational 
analysis shows that cross-talk has an effect on whether the model exhibits bistability. With pathway 
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cross-talk, the behaviour of key proteins switches from transient to sustained activation upon varying the 
duration of the signal. 

We have found only one paper on the application of formal models to cross-talk, ll20lL which contains 
a model of the multiple modes of cross-talk between the EGFR and LIN-12/Notch signalling pathways. 
A discrete, dynamic, state-based model is developed in using the language of Reactive Modules. Model 
checking is used to check the validity of the model and to generate "new biological insights into the 
regulatory network governing the cell fate". However, this work differs from ours because it concerns 
intercellular cross-talk within a multi-cellular model. The application of the term cross-talk to intercellu- 
lar communication is often considered a misnomer amongst the community, and hence we have focussed 
our analysis to intracellular cross-talk in a single-cell model. 



8 Conclusion and Future Work 

In this paper we have explored formal modelling and analysis techniques for signalling pathway cross- 
talk. We note that formal methods, whilst applied frequently in the literature to signalling pathways, 
have been largely ignored in the study of cross-talk. The aim of this paper has been to show that formal 
methods are both well-suited to and a natural choice for this area. 

Our first contribution is a novel categorisation of cross-talk, drawing from examples in the literature. 

The second contribution is the definition of a pathway and pathway composition, based on a set 
of generic pathway modules with internal and external reactions, and renaming and synchronisation 
operations. We can compose pathway independently, or with cross-talk. Furthermore, we also find that 
temporal logic descriptions of behaviours are suitable to detect and characterise cross-talk. The approach 
is illustrated with example pathways. 

The third and final contribution is the application to cross-talk between the TGF-j8/BMP pathway 
and two other pathways. 

Several future directions have been identified. As suggested earlier in the paper, we wish to apply our 
approach to models with a higher level of quantitative detail to make predictions and generate insights 
about the biological effects of cross-talk. Furthermore, we wish to assess how the effect of cross-talks 
differs to that of standard pathway motifs. Hence, is there a reasonable alteration of the pathway model 
(e.g. addition of a feedback loop) that gives the same behaviour as a potential cross-talk? We also 
wish to assess how the effectiveness of pathway intervention techniques (e.g. drugs, gene knockouts) 
changes with the addition of cross-talk. Finally, a larger question is how the temporal ordering of signals 
affects the detectability and behaviour of cross-talk (i.e. do the pathways hold a "biochemical history" 
of signalling events?). 

The PRISM models and CSL properties used in this paper can be found at: 

www . dcs . gla . ac . uk/~radonald/f bt c2010/[ 
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